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The effective baryon-baryon interactions is studied in the refined quark delocalization color screen¬ 
ing model (QDCSM), in which the different quark clusterings are fully taken into account, instead 
of controlling by a variational delocalization parameter e(s) between two 3-quark clusters. The 
symmetry bases are employed to do the calculation, all possible configurations for two quark clus¬ 
ters are considered. The results obtained are very similar to that of QDCSM. It is inferred that 
the delocalization parameter e(s) used in QDCSM is an economic and effective way to describe the 
mixing of quarks between baryons. 


PACS numbers: 12.39.Jh, 13.75.Cs, 14.40.Rt 

I. INTRODUCTION 

Since the advent of quantum chromodynamics (QCD), 
the fundamental theory of strong interaction, it has 
been expected to calculate the properties of baryons 
and baryon-baryon (BB) interactions from QCD directly. 
However, the complexity of the non-perturbative charac¬ 
ter of QCD in the low energy region has hindered this at¬ 
tempt. Although the lattice QCD [jjL Dyson-Schwinger 
equation approach |2j, QCD sum rule Q, chiral perturba¬ 
tion theory |j] and other non-perturbative methods have 
made impressive progresses, the quark model approach 
is still the important method in hadron physics, espe¬ 
cially for BB interaction and multiquark system. The 
quark model describes the experimental data of hadron 
properties and hadron-hadron interactions very well. 

The naive quark model (Glashow-Isgur model 1) suc¬ 
cessfully described the properties of baryons. However 
to extend it to BB interaction, refinements are necessary. 
One approach is to add the Goldstone-boson-exchange 
interaction to the Hamiltonian, it is developed to the 
chiral quark model [6j. After fine-tuning the model pa¬ 
rameters, the model give a satisfactory description of BB 
interactions, where the cr-meson and other scalar mesons 
play an indispensable role. However, the cr-meson, as an 
S -wave resonance of 7T7 t Q, cannot provide the expected 
intermediate-range attraction for nucleon-nucleon inter¬ 
action [§]. In addition, how to realize the chiral partner 
in flavor SU 3 case is still a problem. The outstanding 
feature of the similarity between the nuclear force and 
molecular force can not have an explanation in this me¬ 
son exchange model approach. The another approach 
is the quark delocalization color screening model (QD¬ 
CSM) |9|. Two ingredients, quark delocalization and 
color screening, are introduced to enlarge the model space 
and modify the interactions of quark-pair in different 
quark clusters. This model is based on the same idea 
as Heitler and London’s approach of hydrogen molecular 
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structure. It explains the similarity between molecular 
force and nuclear force naturally. It also gives a good 
description of the BB interaction. Its predictions on d* 
dibaryon (10] had been confirmed by Celsius-WASA ex¬ 
periments [ill ] , the prediction of NQ dibaryon m is also 
supported by lattice QCD calculation [l3j]. In QDCSM, 
the color screening is used to lower the color confine¬ 
ment potential between two baryons and to facilitate the 
quark delocalization. The further studies show that the 
color screening phenomenology is an effective description 
of the hidden-color channels coupling Q The quark 
delocalization is realized by introducing a delocalization 
parameter e(s), which is determined by the dynamics of 
the system, 

M r ) = (<Mr) + e<M r )) /N{e), 
r ) = (<Mr) + e(t> L { r)) /N(e), 

N(e) = a/ 1 + 2ev + e 2 , v = (4> L \4> R ), 

Generally the un-delocalized single particle wavefunc- 
tions take the Gaussian form, 

\L) = (j) L (r) = (71-& 2 ) 4 e _ ^ T<r+s/2) , 

\R)=Mr) = (^ 2 ) _l e-^ (r - s/2)2 (1) 

In fact, the full delocalization can be realized by 
configuration mixing, i.e., taking all the configurations 
\L e ),\L 5 R),\L 4 R 2 ),\L 3 R 3 ),\L 2 R 4 ),\LR 5 ) and |i? 6 } into 
consideration for BB interaction, under the two cluster 
approximation. The mixing of the above seven config¬ 
urations is determined by the dynamics of the system 
directly and no longer by one delocalization parameter e. 
In this way, the model space is larger than that of QD¬ 
CSM. FI. Stancu and L. Wilets had pointed that all the 
configuration are important in their studies of nucleon- 
nucleon interaction |l5| . So it is interesting to study the 
BB interactions with the configuration mixing method 
and to check the validity and efficiency of delocalization 
method in QDCSM. For simplicity, the first version of 
QDCSM, where the interaction between quark pair con¬ 
sists of color confinement and one-gluon-exchange only, 
is employed in the present study. 
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TABLE I: The allowed K for interesting sets of a. The indices 
of the symmetries [v\, [/i], [/] are: l-[6]; 2-[51]; 3-[42]; 4-[33]; 
5-[411]; 6-[321]; 7-[222]; 8-[3111]; 9-[2211], 


a 




K 







201 

144 

264 

324 

344 

394 364 354 




210 

143 

263 

323 

343 

393 

363 

353 




203 

144 

344 









-400 

141 

341 









-1*2 

143 

146 

233 

232 

235 

236 

263 

266 

265 

323 


322 

325 

353 

356 

352 

355 

343 

346 

363 

366 


365 

396 

412 

433 

436 

432 

435 

486 485 

476 


TABLE II: the allowed spatial symmetry for given m, n. 

[i /] L 6 L 5 R L 4 R 2 L 3 R 3 L 2 R 4 LR 5 R 6 
[ 6 ] 1 1 1 1 1 11 

[51] 011 1 110 

[42] 0 0 1 1 10 0 

[33] 0 0 0 1 0 0 0 

TABLE III: The number of symmetry bases for interesting 
physical channels. 

_YU_201 210 203 -400 -1^2 

coupling channels 27 27 10 10 96 


This paper is organized as follows. In section II, the 
model Hamiltonian, the symmetry bases and the calcu¬ 
lation method are described. The results are given in 
section IV and a summary is given in the last section. 


II. MODELS AND BASES 

The model Hamiltonian is the same as that of QDCSM, 
which is described in detail in Ref. [|J. Here we only write 
down the Hamiltonian for 6-quark system, 


H( 6 ) = 


6 

E 

i =1 


2 6 

K + E Vy-T'i 6), 

* i<j= 1 


V■ ■ — V c 4- V c ' 

V IJ — v l0 -|- 


V* = —a r X, ■ A., 


G A, - A ; 

v ij u s 4 


1—e 


1 n xt \ 

— - 2 6{Tij) 
1 11 4 


1 


1 


4 <Ji 


r 2 3 miTO, 


To do the calculation, the symmetry bases are more 
convenient. The symmetry bases are the group-chain 
classification bases, which are defined as follows jH| 




■ Kmn 


(0 = 


[v]L m R n 

[cr }W[/j]P[f]YIJ MiMj 


( 3 ) 


where a = ( YIJ ) are the strong interaction conserved 
quantum numbers: strangeness, isospin and spin. K 
denotes the intermediate quantum numbers, [v], [p], [/]. 
[v \, [/r], [a ], [/] represent the symmetry of orbital, spin- 
flavor SU( 6 ), SU(3) color and flavor, [cr] = [222] is fixed 
due to the color singlet requirement. For some inter¬ 
esting sets of quantum numbers a, the allowed interme¬ 
diate quantum numbers K are listed in Table [H Table 
El presents the allowed spatial symmetries [v\ for given 
m, n. The number of symmetry bases for given a can 
be obtained from the combination of the allowed K and 
m,n, which is given in Table EH 


For illustration, two bases of the ten bases for channel 
a = ( YIJ) = 203 are shown below, 


<T>( 203 ) ( n 6) — 
v 144 5lW i — 

*(203) / 6\ _ 
^344 24 ) ~ 


[6 JL 5 # 1 \ 

[222] 1 [33] 1 [33] 20303 / ’ 

[42] L 2 1? 4 \ 

[222] 1 [33] 1 [33] 20303 / ' 


where Mi,Mj take their maximum values, the eigenen- 
ergies are independent of 

To study the BB interaction, the Schrodinger equation 
(2f)or 6-quark system has to be solved, 


if i,j occur in the 
same baryon orbit, 

if i,j occur in 
different baryon orbits. 


H( 6 )V a {q b ) = E a ^ a (q b ), 


( 4 ) 


where the eigen wavefunction is the linear combina¬ 
tion of ^Kmnil 6 ) under the cluster approximation, 


T“( g 6 ) = £ 


''- 6 


" Kmn*Kmn 


(q b ). 


( 5 ) 


Kmn 


By using Eq.©, Eq.Jf]) becomes 


E [<*2'|tf(6)|*2> - E a (^\n)] <£(9 6 ) = 0, (6) 


where k stands for Kmn. The eigen energy of the system 
can be obtained by solving the generalized eigen equa¬ 
tion. the calculation of the 6-quark Hamiltonian matrix 
elements on the symmetry basis is performed by the well 
known fractional parentage expansion technique fl6| . 

m\H \$£> = E(®X*fcl ai k i,a 2 k 2 )(a[h[,a' 2 k' 2 m,) 

(aiki\a[k[)(a 2 k 2 \H 5 e \a 2 k 2 ). (7) 

where (aikila^k 1 ^ is the four quark overlap. 
(«2 k 2 1 £^ 56 1 ot' 2 k 2 ) is the two body matrix element 
and H$q represents the two-body operator for the last 
pair. |aifci, a 2 k 2 ) and (a^, a 2 k 2 \$%,) are the 

total coefficients of fractional parentage. All the needed 
coefficients can be obtained from the Chen’s book E3- 
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FIG. 1: The effective potentials for NN with YIJ = 201. ’cc’ 
stands for configuration mixing calculation. 

TABLE IV: The parameters used in the calculations. 
m(MeV) m s (MeV) 6(fm) a s a c (MeV-fm -2 ) /x(fm -2 ) 
313 633.76 0.603 1.54 25.13 0 or 1.0 


III. RESULTS 

The effective potentials between two baryons in the 
u, d , s 3-flavor world are calculated in the framework of 
the refined QDCSM. The effective potential is defined as 

V e = E 6 (s) - E 6 (S = oo) (8) 

where Eq(s) is the eigen-energy of 6-quark system with 
separation s. To save space, only the results of several 
interesting states are given below. The needed model 
parameters are all taken from QDCSM, which is listed in 
Table HVl 

Before presenting the results, we discuss the problem 
of numerical calculation first. When the separation be¬ 
tween two clusters is very small, the difference between 
L and R goes to vanishing, the set of bases of the system 
will be over complete. To remove the spurious bases, the 
eigen method is used. First the overlap matrix of the 
system is diagonized, the eigenvectors corresponding to 
zero eigenvalues are dropped. Then the Hamiltonian ma¬ 
trix is reproduced on the remained eigenvectors of over¬ 
lap matrix with non-zero eigenvalues. At last the new 
Hamiltonian matrix is diagonized to get the eigen-energy 
of the system. 

Another problem needs to mention is associated with 
two configurations L 6 and R 6 . Physically, these two con¬ 
figurations are the same, and the energies of these two 
configurations do not change with the separation between 
two clusters. In fact, these two configurations almost do 
not contribute to the effective potentials. 

Figs. 1-5 give the effective potentials between two 
baryons with quantum numbers YIJ = 201, 210, 203, 



FIG. 2: The same as Fig. 1 for NN with YIJ = 210. 

— 1^2, —400. For comparison, the results for // = 0 fm -2 
(without color screening) and the results of QDCSM are 
also shown in these figures. Clearly, the effective poten¬ 
tials go to zero when the separation becomes large (e.g., 
s = 3 fm), due to characteristic of confinement and one- 
gluon-exchange potentials and the Gaussians used as the 
wavefunctions of single-particle. All figures also show the 
differences between the present approach and QCDSM 
gradually disappear when the separation becomes large, 
because the color confinement pushes the energies of the 
hidden color configurations, [&\R 5 L l , [6 ]R 4 L 2 , [6]1? 2 L 4 , 
[6}R 1 L 5 , [42 }R 4 L 2 , [42]i? 2 L 4 , higher. 

Fig. 1 shows the effective potentials between two nu¬ 
cleons with quantum numbers YIJ = 201. If there is 
no color screening effect, i.e., n = 0 fm~ 2 , then the in¬ 
termediate range attraction is missing both in present 
work and in QDCSM, and QDCSM have a larger re¬ 
pulsive core. Taking into account of the color screening 
effect, the intermediate range attraction is obtained in 
both approaches. The attraction in QDCSM is about 
half of that in the present approach. This is reasonable, 
because the variational space in the present approach is 
larger than that in QDCSM. Even though all the config¬ 
urations, \L 6 ),\L 5 R),\L 4 R 2 ),\L 3 R 3 ),\L 2 R 4 ),\LR 5 ) and 
|1? 6 ) are included in QDCSM, the percentages of different 
configurations are controlled by one variational parame¬ 
ter e, which is determined by system dynamics. Whereas 
in the refined QDCSM the dynamics of system fix the 
configuration mixing with 7 coefficients (see Table V). 
By increase the color screening parameter, e.g., fi = 1.6 
fm -2 , QDCSM can produce almost the same effective 
potentials as the present approach. 

The effective potentials for NN with YIJ = 210 are 
shown in Fig. 2. The results of comparison between two 
approaches are similar to that of NN with YIJ = 201. 

For YIJ = 203 A-A channel, attractions appear in 
all cases (see Fig. 3). The color screening introduces 
much stronger attraction. With the same color screening 
parameter, the present approach give a little stronger 
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TABLE V: The relative probabilities of different configurations in the states for YIJ = 203 with /x = 1.0 fm 2 



s = 4 fm 

s = 2.6 fm 

s = 2 fm 

s = 1.1 fm 

s = 0.5 fm 


This work QDCSM 

This work QDCSM 

This work QDCSM 

This work QDCSM 

This work QDCSM 

[6 \R 6 

0.00 

0.00 

0.33 

0.00 

0.27 

0.01 

0.00 

0.02 

0.00 

0.00 

[^L 1 

0.00 

0.00 

0.01 

0.00 

0.47 

0.14 

0.00 

0.19 

0.05 

0.11 

[6 \R 4 L 2 

0.00 

0.00 

0.04 

0.11 

0.74 

0.61 

0.25 

0.67 

0.08 

0.59 

[6 ]R 3 L 3 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

[6 \R 2 L 4 

0.00 

0.00 

0.04 

0.11 

0.74 

0.61 

0.25 

0.67 

0.08 

0.59 

[6 }R 1 L 5 

0.00 

0.00 

0.01 

0.00 

0.47 

0.14 

0.00 

0.19 

0.05 

0.11 

m 6 

0.00 

0.00 

0.33 

0.00 

0.27 

0.01 

0.00 

0.02 

0.00 

0.00 

[42 }R 4 L 2 

0.00 

0.00 

0.02 

0.07 

0.12 

0.09 

0.01 

0.00 

0.01 

0.00 

[42] R 3 L 3 

4.00 

4.00 

3.79 

3.29 

0.61 

0.36 

0.01 

0.00 

0.00 

0.00 

[42]7? 2 L 4 

0.00 

0.00 

0.02 

0.08 

0.12 

0.09 

0.01 

0.00 

0.01 

0.00 


attraction than QDCSM in the short-range part. In¬ 
creasing the color screening parameter in QDCSM, the 
curve will move downward, getting close to the results of 
the present approach. Due to the strong attraction, the 



FIG. 3: The same as Fig. 1 for NN with YIJ = 203. 



FIG. 4: The same as Fig. 1 for NO with YIJ = -1±2. 


bound state with respect to two A’s can be formed. Tak¬ 
ing into account the four-body decay channel, NNtttt , a 
resonance appears in this channel, which was observed 
in WASA-at-COSY experiments (llj. To show the rea¬ 
son for the similarity between two approaches in detail, 
the relative probabilities of different configurations in the 
given states, the ratios of probability of a configuration to 
the configuration [6}R 3 L 3 , are tabulated in Table V. Due 
to the symmetry between left and right Gaussians, the 
probabilities of the configuration [ v]R m L n in the states 
are identical to that of [v\R n L m . From this table, we 
can see that the same physical states are obtained in 
the two approaches when the separation is large enough 
(s = 4 fm). With the decreasing of the separation, the 
differences between the two approaches appear gradu¬ 
ally. However, the main configurations in the states are 
kept the same, [6]I? 3 L 3 , [42]I? 3 L 3 and [6 \R 4 L 2 , [6]I? 2 L 4 . 
Even as the separation becomes very small, the superfi¬ 
cial large difference between two approaches will be re¬ 
duced remarkably because of the large overlap among the 
configurations [6 ]R m L n with different m,n, due to the 
fact that there is no much difference between the right 
and left gaussians. 

The effective potentials of NO with YIJ = —1^2 are 
given in Fig. 4. Again, attractions appear in all cases, 
the color screening increases the attraction about 100 
MeV at the short-range part. A resonance is expected to 
be formed in this channel because of the deep attraction. 
The dynamical calculation of quark model and the lattice 
QCD calculation supported this result. 

Di-Q as a possible dibaryon candidate was predicted 
in 1990 pH and also proposed by Li et al. in 2001 fl9l ]. 
From Fig. 5, we can see that the effective potential be¬ 
tween two fi’s has a mild attraction in QDCSM with 
color screening fj, = 1 fm~ 2 . However, a rather strong at¬ 
traction in the short-range part is obtained in the present 
approach with the same /x. Taking into account the fact 
that the two Q’s can form a weak bound state in QD¬ 
CSM, it is expected that the S2S2 with YIJ = 400 is a 
good dibaryon candidate. 
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FIG. 5: The same as Fig. 1 for QQ. with YIJ = —400. 


IV. SUMMARY 

By enlarging the model space, i.e., taking into account 
of all possible configurations under the constraint of two 


clusters, the effective potentials are obtained for the BB 
systems. In most cases, the potentials obtained in the 
present approach are similar to the ones obtained in QD- 
CSM, except at the short-range part, where lower po¬ 
tentials obtained in this approach. The results show the 
validity of the delocalization in QDCSM. 

From the results for NN channels, the configuration 
mixing themselves does not lead to the intermediate- 
range attraction unless the color screening effect is in¬ 
troduced. The possible dibaryon candidates based on 
the present results are similar to that of the previous 
work jlH ]. Maybe fwill have a little large binding en¬ 
ergy in the present approach. To obtain more reliable 
results for dibaryons, the dynamical calculation is indis¬ 
pensable, which is in progress. 


Acknowledgments 

This work is supported partly by the National Sci¬ 
ence Foundation of China under Contract Nos. 11175088, 
11035006, 11265017 


[1] S. Aoki, PoS (Hadron 2013) 020 (arXive:1402.3059) and 
references therein. 

[2] H. L. L. Roberts, L. Chang, I. C. Cloet and C. D. 
Roberts, Few Body Syst. 51, 1 (2011) and references 
therein. 

[3] L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. 
Rep. 127, 1, (1985); H. X. Chen, A. Hosaka and S. 

L. Zhu, Mod. Phys. Lett. A23, 2234 (2008); H. X. 
Chen, E. L. Cui, W. Chen, T. G. Steele and S. L. Zhu, 
arXiv:1410.0394 

[4] S. Scherer, Adv. Nucl. Phys. 27, 277 (2003); V. Bernard 
and U. Meisner, Ann. Rev. Nucl. Part. Sci. 57, 33 (2007). 

[5] N. Isgur and G. Karl, Phys. Rev. D 18, 4187 (1978); 19, 
2653 (1979); 20, 1191 (1979). 

[6] Y. Fujiwara, C. Nakamoto and Y. Suzuki, Phys. Rev. C 
54, 2180 (1996); Z. Y. Zhang, Y. W. Yu, P. N. Shen et al, 
Nucl. Phys. 625, 59 (1997); A. Valcarce, H. Garcilazo, 
F. Fernandez and P. Gonzalez, Rep. Prog. Phys. 68, 965 
(2005) and references therein. 

[7] M. Ablikim et al. (BES Collaboration), Phys. Lett. B 
598, 149 (2004). 

[8] N. Kaiser, S. Grestendorfer, and W. Weise, Nucl. Phys. 
A 637, 395 (1998); E. Oset, H. Toki, M. Mizobe, and T. 
T. Takahashi, Prog. Theor. Phys. 103, 351 (2000); M. 

M. Kaskulov and H. Clement, Phys. Rev. C 70, 014002 
(2004). 

[9] F. Wang, G. H. Wu, L. J. Teng, and T. Goldman, Phys. 
Rev. Lett. 69, 2901 (1992); G. Wu, L. Teng, J. L. Ping, 
F. Wang, and T. Goldman, Phys. Rev. C 53, 1161 (1996); 
J. L. Ping, F. Wang, and T. Goldman, Nucl. Phys. A657, 
95 (1999); G. H. Wu, J. L. Ping, L. J. Teng, F. Wang, 
and T. Goldman, Nucl. Phys. A673, 279 (2000); H. R. 
Pang, J. L. Ping, F. Wang, and T. Goldman, Phys. Rev. 
C 65, 014003 (2001). 


[10] T. Goldman, K. Maltman, G. J. Stephenson, K. E. 
Schmidt, and F. Wang, Phys. Rev. C 39, 1889 (1989); 
J. L. Ping, F. Wang and T. Goldman, Nucl. Phys. A 
688 , 871 (2001); J. L. Ping, H. R. Pang, F. Wang and T. 
Goldman, Phys. Rev. C 65, 044003 (2002); J. L. Ping, 
H. X. Huang, H. R. Pang, F. Wang, and C. W. Wong, 
Phys. Rev. C 79, 024001 (2009). 

[11] M. Bashkanov et al. (CELSIUS-WASA Collaboration), 
Phys. Rev. Lett. 102, 052301 (2009); P. Adlarson et al. 
(WASA-at-COSY Collaboration), Phys. Rev. Lett. 106, 
242302 (2011); P. Adlarson et al. (WASA-at-COSY Col¬ 
laboration), Phys. Lett. B 721, 229 (2013); P. Adlarson 
et al. (WASA-at-COSY Collaboration), Phys. Rev. Lett. 
112, 202301 (2014). 

[12] T. Goldman, K. Maltman, G. J. Stephenson, K. E. 
Schmidt and F. Wang, Phys. Rev. Lett. 59, 627 (1987); 
H. R. Pang, J. L. Ping, F. Wang, T. Goldman and E. G. 
Zhao, Phys. Rev. C 69, 065207 (2004). 

[13] F. Etminan, H. Nemura, S. Aoki, et al. [HAL QCD Col¬ 
laboration], Nucl. Phys. A 928, 89 (2014). 

[14] H.X. Huang, P. Xu, J. L. Ping and F. Wang, Phys. Rev. 
C 84, 064001 (2011). 

[15] FI. Stancu and L. Wilets, Phys. Rev. C 36, 726 (1987); 
38, 1145 (1988); 40, 1901 (1989). 

[16] F. Wang, J. L. Ping and T. Goldman, Phys. Rev. C 51, 
1648 (1995); 3346 (1995). 

[17] J.Q. Chen et al, Tables of the Clebsh-Gordan, Racah and 
Subduction Coefficients of SU ( n) Groups, 1987 (World 
Scientific, Singapore); Tables of the SU(mn ) D SU(m) x 
SU(n) Coefficients of Fractional Parentage, 1991 (World 
Scientific, Singapore). 

[18] V. B. Kopeliovich, B. Schwesinger, and B. E. Stern, Phys. 
Lett. B 242, 145 (1990). 

[19] Q. B. Li, P. N. Shen, Z. Y. Zhang, and Y. W. Yu, Nucl. 










6 


Phys. A 683, 487 (2001) 



